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ABSTRACT 

One of the most compelling tasks of modem cosmology is to constrain the expansion his¬ 
tory of the Universe, since this measurement can give insights on the nature of dark energy 
and help to estimate cosmological parameters. In this letter are presented two new measure¬ 
ments of the Hubble parameter H(z) obtained with the cosmic chronometer method up to 
z ~ 2. Taking advantage of near-infrared spectroscopy of the few very massive and passive 
galaxies observed at z > 1.4 available in literature, the differential evolution of this popula¬ 
tion is estimated and calibrated with different stellar population synthesis models to constrain 
H(z), including in the final error budget all possible sources of systematic uncertainties (star 
formation history, stellar metallicity, model dependencies). This analysis is able to extend sig¬ 
nificantly the redshift range coverage with respect to present-day constraints, crossing for the 
first time the limit at z ~ 1.75. The new H(z) data are used to estimate the gain in accuracy on 
cosmological parameters with respect to previous measurements in two cosmological models, 
finding a small but detectable improvement (~5 %) in particular on and wq. Finally, a 
simulation of a Euclid-like survey has been performed to forecast the expected improvement 
with future data. The provided constraints have been obtained just with the cosmic chronome¬ 
ters approach, without any additional data, and the results show the high potentiality of this 
method to constrain the expansion history of the Universe at these redshifts. 
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1 INTRODUCTION 

Since the pion eering works of iRiess et all d 1998h and 
IPerlmutter et alJ i ll 9991) . the comprehension of the acceler¬ 
ated expansion of the Universe has become a key question in 
cosmology. Amongst many proposed approache s to tackle this 
open issue (for a comprehensive review, see IWeinbere et ali 
|2013h . the most common ones rely on probes which are able to 
disentangle the evolution with redshift due to the expansion of 
the Universe, and the intrinsic evolution of the considered probe; 
for this purpose, generally probes with “standard” properties have 
been considered, such as standard candles (Supernovae Type la, 
SNe) and standard rulers (Baryon Acoustic Oscillations, BAO). In 
this context, each analysis has its own strength and weakness, but 
it is now common understanding that it is of primary importance to 
explore various possibility at the same time, to keep under control 
the systematic effects biasing each single probe, and to cut down 
the errors on cosmological parameters by combining the different 
probes. 

As suggested by llimenez & Loebl ( 2002 1. an alternative 
method to constrain the expansion history of the Universe is based 
on the study of the redshift evolution of “cosmic chronometers”. In 
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this approach, the relative age of passive galaxies dz/dt can be used 
as standard clocks , since H(z) = —1/(1 +z)dz/dt. The key issues 
of this method are the identification of a proper population of galax¬ 
ies to trace homogeneously the differential age evolution of the 
Universe, and the reliance on theoretical models to estimate galaxy 
ages, and hence the degeneracy between parameters. As poin ted out 
by many authors (for an extensive review, see lRenzin i OOOfil ). mas¬ 
sive and passive early-type galaxies represent the ideal population 
to address the first point, having assembled their mass at high red¬ 
shift and being passively evolving since then. To solv e the second 
point, in previous papers dMoresco et al.ll20Til . l20 1 2r3) an improve¬ 
ment of the technique was proposed based on the study of a direct 
observable of galaxy spectra, the 4000 A break (D4000). which 
is known to be strictly correlated with galaxy age and metallicity, 
and less significantly dependent on star formation history (SFH); 
since the assumption of a linear relation between £>4000 and age 
of a galaxy has been proven to b e an extremely good approxima¬ 
tion for most £>4000 regimes (see lMoresco et al.ll2.01 ill 2012/1) . i.e. 
£>4000 =A(SFH,Z/Zq) • age + 5, it is possible to redefine H(z) as: 

—sk A{SFHZIZ °)jMm "> 

The strength of this new equation is that now the dependence on 
statistical and systematic effects has been decoupled: the factor 
dz/dD4000 is only dependent on observables, while all the depen- 
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dence on degeneracy between parameters or assumption of stel¬ 
lar population synthesis (SPS) models is contained in the factor 
A(SFH ,Z/Zq), which can be calibrated on SPS models. Its draw¬ 
back is that it relies on the estimate of the stellar metallicity of the 
population, or at least on having a reasonable proxy of it, to be used 
as a prior. All the possible systematic error deriving from these as¬ 
sumptions (both on models and on the galaxy properties) have been 
considered and taken in consideration in the total error budget, as 
discussed in the following section. In iMoresco et al.l d2012al) was 
demonstrated that this new approach relying on the spectroscopic 
differential evolution of cosmic chronometers is robust against 
the choice of SPS models, and not significantly dependent on the 
SFH assumptions (also because the selection criteria leaves small 
space for prolonged SFHs). In that paper, eight new H(z) measure¬ 
ments have been provided in the redshift range 0.15 < z < 1.1, 
significantly extending the redshift coverage and precision o f pre¬ 
vious similar analysis dSimon et al.l 120051 : Istem et alJl2010h : the 
potential of this new method in comparison with more standard 
probe s has been studied by many authors (e.g . see Jimgnezetalj 


probes nas been studied by many authors (e.g. see Jimenez et at. 
20121: IMoresco et alJl2012bl: Izhao et all 1201 1: IWang et alj|2012l: 


Rienier-Sdrensen et al.l 20131 ; Verde et alJ 20141) . which demon¬ 


strated how the cosmic chronometers method can be competitive 
for many aspects with SNe and BAO in constraining cosmological 
parameters. 

In this letter two new H(z) measurements up to z ~ 2 are re¬ 
ported, taking advantage of various literature data of high redshift 
(z > 1.4) massive and passive galaxies with near-infrared spec¬ 
troscopy covering the restframe range around the D4000. These 
data, never used before for this scope, represent a golden mine 
to extend present-day estimates of H(z) to a redshift range never 
approached before with this technique, and matched only by the 
measurement of the BAO in the Lya forest of BOSS quasars by 
iBusca et~al1 d201 3l) . Throughout this letter the new data are pre¬ 
sented, and the gain in accuracy for different cosmological parame¬ 
ters discussed, providing also forecasts of improvement which may 
be reached with future surveys. 


2 DATA 

Different literature data providing spectroscopical analysis of high 
redshift (z > 1.4) massive and passive galaxies are considered. 
The most massive, passive and with the shortest SFH sample has 
been selected from each analysis: these galaxies represent the most 
appropriate probe to sample the expansion history of the Uni¬ 
verse, mapping the “eldest crust” population at each redshift in¬ 
dependently of the selection criteria. To apply the method, it has 
been considered the definition of the break with narrower bands 
(D4000„,3850-3950 and 4000-4100 A) introduced bv lBalogh et al.l 
J 1999t) , being less sensitive to potential reddening effects. While 
some works provide such measurement, for the others the spectra 
have been analyzed, obtaining for the first time D4000 n estimates 
for these samples. 

In the following, the samples used, and how they have been se¬ 
lected, are briefly presented. For the full description these analyses, 
one can refer to the original papers. 

• ISaracco et alJ d2005l) . This work presents seven massive 
evolved galaxies spectroscopically identified in the MUNICS sur¬ 
vey; the spectral analysis confirmed their passive nature, with short 
SFH ( x < 0 .3 Gyr) and high masses (M/Mq > 10 11 ). 

• iKriek et all 1 20091) . They analyze the spectrum of a quies¬ 
cent, ultra-dense galaxy at z=2.1865 obtained with the Gemini 
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Figure 1 . Redshift (left plot) and stellar mass (right plot) histograms for the 
high-z passive galaxies analyzed. The red dotted lines show how the sample 
has been divided into three redshift bins. 


Near-Infrared Spectrograph. The spectral fitting provided a low 
star formation rate (present SFR<1% of past SF), and large mass 
(M/ Mq, > IQ"), _ 

• lOnodera et al .1 d20k3) . In this work. 18 passive BzK galaxies in 
the COSMOS field at z > 1.4 were spectroscopically observed and 
analyzed; the fit to their SED provided stellar masses M/Mq = 4 — 
-40 x 10 10 . The most passive sample is further selected, discarding 
galaxies with 24 /rm detection, prolonged SFH (t >= 0.3 Gyr) and 
high HS values (EW(H5 ) <-3A). 

• iGobat et alJ j 20131) . From HST/WFC3 slitless spectroscopic 
observations of the cluster Cl 11449+0856, 6 quiescent members 
were identified at z = 2; four of them have a high enough S/N to al¬ 
low the measurement of D4000, and their stellar masses were mea¬ 
sured in StrazzulloetalJj 20131) (M/Mq > 1.2 x 10 11 ). 

• iKrogager et al.1 j2013l) . In this work, they presented deep 
HST/WFC3 grism spectroscopy for a sample of 14 galaxies at z ~ 2 
in the COSMOS field. Among them, 7 are quiescent galaxies with 
measured D4000 and masses M/Mq > 0.6 x 10 11 . Due to the low 
resolution and poor sampling of the grism data, the D4000 has been 
measured in its definition with larger bands (D4000,,., 3850-3950 
and 4000-4100 A. lBruzualll983l) . and to compare it with the other 
data SPS models are used to obtain the required conversiorj]]. 


In this way, 29 galaxies are selected, which represent the most mas¬ 
sive, red and passive galaxies found so far at these redshifts. On 
average, the formation redshifts (estimated from the best-fit to both 
galaxy spectra and photometry) are 3.5 < Zf < 5, with a spread in 
formation age < 0.45 Gyr, which is much smaller than the asso¬ 
ciated error, therefore avoiding possible progenitor bias issues. In 
Fig.m are shown the redshift and stellar mass distributions of the 
final sample. The scatter in D4000,, between the analyzed galax¬ 
ies (0/54000 ~0.2-0.25) is statistically consistent with the properties 
of the considered population, with the D4000,, errors and the es¬ 
timated SFH contributing in similar fractions to it. To apply the 
method, the sample has been further divided into three redshift 
bins to approximatively preserve the same number of galaxies in 
each bin, as shown in the figure. This choice is a trade-off aimed to 
maximize both the number of H(z) measurements and to minimize 


1 Note that the conversion factor from D4000,,. to D4000,, in this range of 
values is ~ 12%, and does not significantly depend on the assumed model 
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Figure 2. Upper plot: D4000,, -z relation for th e galax ies considered in this analysis (black points), their weighted average (red points) and the median relation 
obtained at lower redshifts in iMoresco et alJ12012 a l in the same mass regime. Intermediate plot: Hubble parameter measurements obtained (red points), 
compared with literature data. Filled points have been obtained with D4000-based cosmic chronometers, while the open points from the relative age based 
method. The theoretical relation is not a best-fit to the data, but the 11(7.) obtained assuming a Planck fiducial cosmology (Cl m = 0.315, Ho=67.3 km/sec/Mpc, 
Ipianck Collaborationll2014h . The best fits to the data for a open ACDM cosmology and for a flat wCDM cosmology are shown respectively with a dotted and 
dashed line. Lower plot: Simulated H(z) data from a Euclid-like survey. The inner errorbar represent the statistical error, while the outer enorbar the total error. 


the relative errorbars. The Z94000,, data have been then averaged 
in the previously discussed redshift bins using a weighted mean, 
to take into account also the error associated to each sample, since 
some points have a much smaller error than others; the results are 
shown in the upper plot of Fig. [2] In the figure it is also shown 
the median D4000 n -z relation obtained in iMoresco et alj d2012al) 
in the same mass range (10* 1 < M/Mq < 10 1125 ), to give an in¬ 
sight of the overall evolution between z ~ 0 to z ~ 2. The last point 
of lMoresco et al.1 d2012ah will be also used as connection point be¬ 
tween the low-redshift and the high-redshift analysis, to estimate 
the first H(z) point. 


3 METHOD AND ANALYSIS 

The Hubble parameter H(z) is evaluated following eq.Q] where the 
relevant quantities are dz/dD4000 n and A(SFH,Z/Zq). To obtain 
the first one, it has been estimated the evolution between the last 
D4000,, point provided in lMoresco et alj d2012ah and the first point 
of this analysis, and between the second and third points of this 
analysis. In this way, it was possible to provide two new H(z) mea¬ 
surements. The parameter A(SFH,Z/Zq) has been estimated from 
stellar p opulation synthesis models, as described in lMoresco et al.l 
d2012al) . Briefly, given the properties of the selected population, are 
considered models with an exponentially delayed SFH with differ¬ 
ent x (up to x=0.3 Gyr, characteristic of the selected sample) and 
a stellar metallicity in the range 0.75 < Z/Zq < 1.25. This last 
choice is well justified both on observational and theoretical ba¬ 
sis. On one hand, for the selected sample only broad constraints 
on metallicity are provided, with a general agreement on a near 


Table 1. Hubble parameter measurements. 


z 

H(z) [km/sec/Mpc] 

1.363 

160 ±33.6 

1.965 

186.5±50.4 


solar (or slightly subsolar) metallicity. The most precise constraint 
so far for this pop ulatio n of very mas sive and passive galaxies has 
been obtained by lOnodera et al.l (I2014I) . finding from the stacked 
spectrum of 24 galaxies a metallicity in the range 0.024± 20%. 
On the other hand, these passive systems have already exhausted 
their gas reservoir, and therefore a negligible evolution is expected 
with respect to their low-redshift counterparts, which show an al¬ 
most solar metallicity. Therefore, a conservative range of ±25% 
around solar metallicity is a ssumed. Two ste llar population synthe- 


lar po 

sis models have been used: Bruzual & Charloj d2003l) (BC03) and 
iMaraston & Strombackl d201 lh (MaStro). They encompass differ¬ 
ent assumptions for stellar evolutionary models, treatment of the 
asymptotic giant branch phase, and method used to compute the 
integrated spectra, and so can be used to quantify the dependence 
of the results on the models. The parameter A(SFH,Z/Zq) is es¬ 
timated in both models, and the difference in the resulting H(z) 
values has been summed in quadrature to the total error. 

In Tab.|T|and Fig.[2]are reported the two new H(z) constraints, 
and for comparison are shown also the values available in liter¬ 
ature obtained from “cosmic chronometers” approach, both from 
relative ages or D4000. As it is possible to see, this analysis for 
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Table 2. Marginalized cosmological parameters constraints obtained with 
cosmic chronometers in the two considered cosmologies. The errors re¬ 
ported correspond to I -a confidence levels. 


oACDM 

O.M 

Cl A 

Q. k 

D4000 H(z) 

0-4*±P? 

0 93+ 0 07 

—0.48 

—0 43+ 0 - 87 
u -^ 3 —0.23 

D4000 H(z) + new data 

o-4i±8:| 2 

0 74+ 0 - 26 
u - / ^-0.34 

—0 15+ 0 - 68 

U ' 1 —0.45 

all H(z) 

0 27 +0 ' 19 
u - z/ — 0.19 

0 55 +0 ' 27 
—0.25 

-0 19+ 0 - 44 
uiy -0.46 

all H(z) + new data 

0 26 +0-18 
u - zo -0.18 

0 54+0- 26 
u - J ^-0.24 

0 21 +0-42 
u ’ Zi -0.45 

all H(z) + new + sim data 

0 ?4 +0-13 
—0.13 

0 51 +0.22 
u - :)1 -0.21 

0 25 +0-34 
u - zd -0.35 

fwCDM 

Ho [km/sec/Mpc] 

O.M 

wo 

D4000 H(z) 

67.8^ 2 5 8 

0 32 +0,11 
u - 3Z -0.08 

< -0.78 

D4000 H(z) + new data 

67.3l‘ 2 3 5 

0 31 + 0 - 1 

u-- ,1 _0.08 

< -0.5 

all H(z) 

67+1:1 

n 0+0.67 
—0.06 

1 +0.58 
—0.92 

all H(z) + new data 

66.8 tH 

r\ O+0.06 
u ' —0.06 

-0.96+g|t 

all H(z) + new + sim data 

66t 6 4 3 

0 29 +0 ' 06 
u,z —0.05 

-0.82+>£ 


the first time breaks the current limit at z ~ 1.75, proving a mea¬ 
surement at z ~ 2. These results have been obtained with just a few 
tens galaxies, and shows the potentiality of this approach to provide 
valuable constraints in the range 1.5 < z < 2; in this range, forth- 
coming surveys suc h as Euclid lLaureiis et al.ll20ill) and WFIRST 
(Sperg el et alj2013h . will greatly improve the present-day statistic, 
increasing the number of passive galaxie s at th ese redshifts by at 
least an order of magnitude (see lLaureiis et alj|201 ill . In estimat¬ 
ing the statistical strength of these new data, a simulation of the 
possible new constraints that can be achieved with new data was 
also performed. It has been considered a Euclid-like survey ob¬ 
serving approximatively one thousand spectra of passive galaxies 
in five redshift bins in the range 1.5 < z < 2; the expected error 
is estimated taking into account, to be conservative, only the im¬ 
provement due to the increased statistic. In this way, the systematic 
errors are expected to clearly dominate over the statistical errors 
(see Fig. [2}. The data have been simulated on a Planck fiducial cos¬ 
mology, shifting the simulated points from the theoretical relation 
with a gaussian distribution with dispersion equal to the estimated 
error. These simulated data are also shown in Fig. [2] 

To explore the constraining power of the new H(z) data, I con¬ 
sidered two models, an open ACDM model: 

H(z) = Hq\J £2jvf(l +z) 3 +£ijt(l +z ) 2 +£2a (2) 

where £2 m + Cl A + £\ = 1 and a flat wCDM model: 

H{z) = H 0 ^n. M (l+z) 3 +a A (l+z) 3 ( 1+w °'> (3) 

A standard % 2 analysis is performed as a function of three cos- 
mological parameters, i.e. jflivi. Qa , flu | (Ho=67.3 km/sec/Mpc, 
Planck Co 11 a bo r atior] 1 201 4) and {£2 m,Ho, w o} respectively. The 
cosmological constraints obtained are shown in Tab. [2] and Figs. 
[3] and [4j the values provided have been obtained by marginalizing 
each time the full 3D likelihood on the uninteresting parameters. 
Both the constraints that can be achieved with only D4000 data and 
with all the H(z) measurements obtained from cosmic chronome¬ 
ters are considered, and in both case the results obtained with and 
without the two new H(z) points of this analysis are compared. For 
many parameters, it is found a measurable, even if small, improve¬ 
ment adding the new H(z) measurements, especially for and wo 
(~5%, see also Fig. 0: this is in particular remarkable considering 


oACDM model oACDM model 



H n [km sec~' Mpc £2 m 


Figure 3. Marginalized contour plots for an open ACDM cosmology (four 
upper panels) and in a flat wCDM cosmology (four lower panels). In each 
cosmology are shown 1-0, 2-0, and 3-0 confidence levels (light, medium 
and dark gray respectively) considering only the D4000 data or all the data, 
as highlighted in the captions. Solid lines show the constraints obtained in¬ 
cluding the new H(z) data, while dashed lines show the constraints without 
them; the relative gain in accuracy obtainable with these new data can be 
inferred by comparing the two lines. 


that this analysis has been performed on a very small set of galax¬ 
ies. The forecasts for a Euclid-like survey are more optimistic, with 
a percentage improvement on the error of ~30%, 20% and 15% re¬ 
spectively for O.M, &k an d Cl A in a oACDM cosmology, and ~10%, 
15% and 20% respectively for Q.m, Hq and wo in a flat wCDM cos¬ 
mology. As expected the improvement on the low-redshift parame¬ 
ter Ho is less significant, since the new data primarily constrain the 
evolution of H(z) at high redshift. 


4 CONCLUSIONS 

This analysis shows the potentiality of using the “cosmic 
chronometers” approach to constrain the expansion history of the 
Universe up to z ~ 2. Using spectroscopic observation of a few 
very massive and passive galaxies in the range 1.4 < z < 2.2, I 
provide two new measurements of H(z) at z = 1.36 and z = 1.97, 
improving the current limit at z ~ 1.75. The new data are shown to 
be important to further improve the accuracy especially Q.m and wo 
with respect to current measurements obtained with this technique. 
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Figure 4. One-dimensional marginalized x 2 for cosmological parameters in an open ACDM cosmology (upper panels) and in a flat wCDM cosmology (lower 
panels). In each case it is shown for comparison the constraints obtained with D4000 data (with and without the new lHz) measurements, grey and red dashed 
lines), with all the data (with and without the new H(z) measurements, black and red solid lines), and also with the simulated data added (blue dashed lines). 


The forecasts obtained on a Euclid-like survey confirmed an even 
higher improvement in the accuracy reachable on most cosmo¬ 
logical parameters. It is of particular interest that the provided 
constraints have been obtained with a single cosmological probe. 
Much more precise estimates will be available once different 
probes are combined, as shown e.g. in lMoresco et alj | |2012bi) . 
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